Phonon  counting  and  intensity  interferometry  of  a  nanomechanical  resonator 
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Using  an  optical  probe  along  with  single  photon  detection  we  have  performed  effective  phonon 
counting  measurements  of  the  acoustic  emission  and  absorption  processes  in  a  nanomechanical  res¬ 
onator.  Applying  these  measurements  in  a  Hanbury  Brown  and  Twiss  set-up,  phonon  correlations  of 
the  nanomechanical  resonator  are  explored  from  below  to  above  threshold  of  a  parametric  instability 
leading  to  self-oscillation  of  the  resonator.  Discussion  of  the  results  in  terms  of  a  “phonon  laser”, 
and  analysis  of  the  sensitivity  of  the  phonon  counting  technique  are  presented. 


In  optics,  the  ability  to  measure  individual  quanta  of 
light  (photons)  enables  a  great  many  applications,  rang¬ 
ing  from  dynamic  imaging  within  living  organisms  [1] ,  to 
secure  quantum  communication  [2  .  Landmark  photon 
counting  experiments,  such  as  the  intensity  interferome¬ 
try  performed  by  Hanbury  Brown  and  Twiss  (HBT)  [3] 
to  measure  the  angular  width  of  visible  stars,  have  played 
a  critical  role  in  our  understanding  of  the  full  quantum 
nature  of  light  [4] .  In  this  work  we  use  an  optical  probe 
and  single  photon  detection  to  study  the  acoustic  emis¬ 
sion  and  absorption  processes  in  a  silicon  nanomechanical 
resonator,  and  perform  a  similar  HBT  measurement  to 
measure  correlations  in  the  emitted  phonons  as  the  res¬ 
onator  undergoes  a  parametric  instability  formally  equiv¬ 
alent  to  that  of  a  laser  [5] .  Owing  to  the  cavity-enhanced 
coupling  of  light  with  mechanics,  this  effective  phonon 
counting  technique  has  a  noise  equivalent  phonon  sen¬ 
sitivity  of  tt-nep  =  0.89  zb  0.05  at  an  intracavity  pho¬ 
ton  number  of  order  unity.  With  straightforward  im¬ 
provements  to  this  method,  a  variety  of  quantum  state 
engineering  tasks  employing  mesoscopic  mechanical  res¬ 
onators  will  be  enabled  [6] ,  including  the  generation  and 
heralding  of  single  phonon  Fock  states  [F  and  the  quan¬ 
tum  entanglement  of  remote  mechanical  elements  m- 
Mechanical  resonators  are  currently  employed  in  myr¬ 
iad  technologies,  from  precision  navigation  systems  m 
to  personal  electronics  m-  As  with  objects  at  the 
atomic  scale,  the  properties  of  macroscopic  mechanical 
resonators  are  governed  by  the  laws  of  quantum  mechan¬ 
ics,  providing  fundamental  quantum  limits  to  the  sen¬ 
sitivity  of  mechanical  sensors  and  transducers.  Current 
research  in  the  area  of  cavity-optomechanics  seeks  to  ex¬ 
plore  the  quantum  mechanical  properties  of  mechanical 
systems  ranging  in  size  from  kilogram  mass  mirrors  to 
nanoscale  membranes  m,  and  to  develop  technologies 
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for  precision  sensing  [13]  and  quantum  information  pro¬ 
cessing  nmg.  Measurement  of  the  properties  of  me¬ 
chanical  systems  in  the  quantum  regime  typically  involve 
heterodyne  detection  of  a  coupled  optical  or  electrical 
field,  yielding  a  continuous  signal  proportional  to  the 
mechanical  displacement  amplitude  m-  An  alternative 
method,  particularly  suited  to  optical  read-out,  is  to  uti¬ 
lize  photon  counting  as  a  means  to  probe  the  quantum 
dynamics  of  the  coupled  opto-mechanical  system  [201 12T]. 

Photon  counting  can  be  readily  adapted  to  study  inten¬ 
sity  correlations  in  an  optical  field,  and  has  been  used  not 
only  in  the  astronomical  HBT  studies  of  thermal  light, 
but  also  in  early  studies  of  the  photon  statistics  of  laser 
light  and  single-atom  fluorescence  [4j  [22] .  In  the  field 
of  photon-correlation  spectroscopy,  such  intensity  inter¬ 
ferometry  techniques  have  found  widespread  application 
in  the  measurement  of  particle  and  molecular  motion  in 
materials  [23] .  More  recently,  photon  counting  of  Raman 
scattering  events  in  diamond  has  been  used  to  herald 
and  verify  the  quantum  entanglement  of  a  THz  phonon 
shared  between  two  separate  bulk  diamond  crystals  [9]. 
In  the  case  of  engineered  cavity-optomechanical  systems, 
much  longer  phonon  coherence  times  are  attainable,  al¬ 
beit  at  lower  mechanical  frequencies  (MHz-GHz)  which 
limit  the  temperature  of  operation  and  the  optical  power 
handling  capability  of  such  structures.  Quantum  opti¬ 
cal  schemes  for  manipulation  of  the  quantum  state  of 
motion  in  cavity-optomechanical  systems  thus  rely  on  a 
large  per-phonon  scattering  rate  and  efficient  detection 
of  scattering  events.  In  this  work  we  embed  a  high-Q, 
GHz-frequency,  mechanical  resonator  inside  an  optical 
nanocavity,  greatly  enhancing  the  phonon-photon  cou¬ 
pling  rate  and  channeling  optical  scattering  into  a  pre¬ 
ferred  optical  mode  for  collection,  both  of  which  combine 
to  yield  a  sub-phonon-level  counting  sensitivity  of  the  in¬ 
tracavity  mechanical  resonator  occupancy. 

A  conceptual  schematic  of  the  phonon  counting  experi¬ 
ment  presented  here  is  shown  in  Fig.  Ek  The  device  stud¬ 
ied  in  this  work  consists  of  a  patterned  silicon  nanobeam 
which  forms  an  optomechanical  crystal  (OMC)  [24j  [25] 
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FIG.  1.  a,  Schematic  of  the  phonon  counting  measurement. 
The  finite  element  method  simulations  depict  the  displace¬ 
ment  held  of  the  acoustic  resonance  (top)  and  the  electric 
held  of  the  optical  resonance.  Pump  light  at  optical  detun¬ 
ing  A  =  —  ujm  (A  =  ct;m)  is  indicated  by  the  blue  (red)  arrows 
and  corresponds  to  measurement  of  phonon  emission  (absorp¬ 
tion),  while  optically  resonant  light  is  represented  by  black 
arrows.  The  optical  rehection  is  hltered  to  reject  the  pump, 
then  detected  in  a  Hanbury  Brown  and  Twiss  (HBT)  set-up 
using  two  superconducting  single  photon  detectors  (SPDs), 
with  outputs  used  as  start/stop  pulses  in  a  time-correlated 
single  photon  counting  module.  Left  side  of  the  schematic 
shows  the  resulting  anti-normal-ordered  phonon  correlation 
function  for  blue  detuned  pumping  (A  =  —  cjm),  while  right 
side  of  the  schematic  shows  red  detuned  pumping  (A  =  cjm) 
which  yields  a  normal-ordered  correlation  function,  b,  Mea¬ 
sured  power  spectral  density  S(u)  of  the  acoustic  resonance, 
c,  Normalized  optical  cavity  rehection  spectrum  R.  Pump  de¬ 
tunings  of  A  =  ±ct;m/(27r)  &  ±5.6  GHz  are  indicated  by  the 
red  and  blue  dashed  lines,  respectively,  d,  Transmission  spec¬ 
trum  of  the  hrst  (purple)  and  second  (orange)  optical  hlter, 
with  total  hlter  transmission  plotted  in  black. 


able  to  co-localize  acoustic  and  optical  resonances  at  fre¬ 
quencies  u;m  and  cjc,  respectively.  Finite-element  method 
simulations  of  the  acoustic  and  optical  resonances  are 
shown  at  the  top  of  Fig.  Ek  The  Hamiltonian  describing 
the  interaction  between  the  acoustic  and  optical  modes 
is  given  by  H[nt  =  hgod^a(b  ±  £F),  where  a  (b)  is  the  an¬ 
nihilation  operator  for  the  optical  (acoustic)  mode,  and 
go  is  the  optomechanical  coupling  rate,  physically  repre¬ 
senting  the  optical  frequency  shift  due  to  the  zero-point 
motion  of  the  acoustic  resonator.  This  interaction  mod¬ 
ulates  an  optical  cavity  drive  at  frequency  uj\  to  produce 
sidebands  at  frequencies  u\  ±  co>m,  analogous  to  the  anti- 
Stokes  and  Stokes  sidebands  in  Raman  scattering  and 
corresponding  to  phonon  absorption  or  emission,  respec¬ 
tively.  For  a  system  in  the  resolved  sideband  limit,  where 
c k  (k  is  the  linewidth  of  the  optical  resonance),  the 
density  of  states  of  the  optical  cavity  can  be  used  to  res¬ 
onantly  enhance  either  scattering  process  for  an  appro¬ 
priately  detuned  pump.  In  particular,  applying  a  large 


coherent  pump  red  (blue)  detuned  from  the  optical  cav¬ 
ity  resonance  by  A  =  ujc  —  uo\  =  (A  =  —  cjm)  re¬ 
sults  in  an  effective  interaction  Hamiltonian  of  the  form 
Hint  =  TiG(a)b  ±  a&)  (Hint  =  hG(ab  +  d^)),  where 
G  =  goy/rG  is  the  parametrically  enhanced  optomechan¬ 
ical  coupling  rate  (nc  is  the  intracavity  photon  number 
at  frequency  u)\  due  to  the  pump  laser).  In  this  case,  the 
output  field  annihilation  operator  dout  can  be  shown  to 
consist  of  a  coherent  component  at  frequency  u\  as  well 
as  a  component  at  frequency  c jc  which  is  proportional 
to  b  (fet)  [26]  (see  App.  [Fjfor  details).  Sending  the  cav¬ 
ity  output  through  a  series  of  narrowband  optical  filters 
centered  on  the  cavity  resonance,  as  shown  in  Fig.  Ek 
suppresses  the  pump  so  that  photon  counting  events  will 
correspond  directly  to  counting  phonon  absorption  (emis¬ 
sion)  events  [27] .  Subsequently  directing  the  filter  output 
to  a  Hanbury  Brown  and  Twiss  (HBT)  set-up  in  order 
to  measure  the  second-order  photon  correlation  function 
g(2\r)  [3]  will  then  result  in  a  direct  measurement  of  the 
normally  (anti- normally)  ordered  second-order  phonon 
correlation  function. 

As  described  in  Ref.  m,  the  nanobeam  is  patterned 
in  such  a  way  as  to  support  a  “breathing”  acoustic  res¬ 
onance  at  cjm/27r  =  5.6  GHz  as  well  as  a  fundamen¬ 
tal  optical  resonance  at  a  free-space  wavelength  near 
1550  nm,  with  a  theoretical  vacuum  coupling  rate  of 
go/27r  =  860  kHz.  As  detailed  in  Ref.  [28],  coupling 
to  the  optical  resonance  is  accomplished  via  a  single¬ 
sided  end-fire  coupling  scheme  using  a  lensed  optical 
fiber.  All  device  measurements  presented  here  are  pe- 
formed  at  room  temperature  and  pressure.  The  thermal 
Brownian  motion  of  the  acoustic  resonance  can  be  ob¬ 
served  as  a  Lorentzian  response  centered  around  u)m  in 
the  noise  power  spectral  density  (NPSD)  S(u)  of  the  pho¬ 
tocurrent  of  the  reflected  optical  signal  from  the  cavity, 
as  shown  in  Fig.  [TJd.  The  linewidth  of  this  Lorentzian  is 
given  by  7  =  71  ±  7om,  where  71  is  the  intrinsic  acous¬ 
tic  energy  damping  rate  and  70M  =  ±4 G2  / k  is  the  op- 
tomechanically  induced  damping  rate  due  to  dynamical 
backaction  when  pumping  on  the  red  or  blue  sideband, 
respectively  [26 [.  By  measuring  this  linewidth  as  a  func¬ 
tion  of  nc  for  both  red  and  blue  detuning,  we  can  extract 
7i/27r  =  3  MHz  and  go/Zir  =  645  kHz.  The  optical  reflec¬ 
tion  spectrum  of  the  cavity  is  shown  in  Fig.  [1J3,  with  the 
frequencies  used  for  red  and  blue  detuned  pumping  indi¬ 
cated  by  dashed  lines.  For  this  device,  we  observe  a  total 
optical  energy  decay  rate  of  k/27t  =  817  MHz  and  a  decay 
rate  into  the  detection  channel  of  Ke/ 2tt  =  425  MHz. 

The  sensitivity  of  the  phonon  counting  scheme  pre¬ 
sented  here  depends  upon  the  above-measured  optome¬ 
chanical  cavity  parameters,  the  dark  count  rate  of  the 
photon  detectors,  Tdark,  and  the  residual  pump  laser 
light  which  is  transmitted  through  the  filters.  In  this 
work  we  use  a  cascaded  pair  of  tunable  Fabry-Perot  fil¬ 
ters  with  bandwidth  of  50  MHz  and  free  spectral  range 
of  20  GHz  (see  Fig.  EH  The  single  photon  detectors 
are  fiber-coupled,  WSi-based  superconducting  nanowire 
single-photon  detectors  (SPDs)  [29].  Each  detector  has 
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FIG.  2.  a,  Noise  equivalent  phonon  number  ttnep  versus  in¬ 
tracavity  photon  number  nc  calculated  using  the  measured 
signal  and  noise  count  rates  for  our  current  setup  (black  cir¬ 
cles).  Solid  lines  indicate  the  theoretically  expected  contribu¬ 
tions  due  to  dark  counts  (red)  and  pump  bleed-through  (blue), 
based  on  the  measured  system  efficiency  and  pump  suppres¬ 
sion.  The  blue  dashed  line  indicates  the  expected  pump  bleed- 
through  contribution  with  an  additional  filter,  b,  71nep  versus 
filter-pump  detuning  Amter  for  nc  ~  65,  with  (red)  and  with¬ 
out  (black)  an  additional  C-band  band-pass  filter  inserted. 
The  vertical  green  line  indicates  the  detuning  corresponding 
to  the  data  from  a,  and  the  horizontal  black  line  indicates  the 
expected  limiting  sensitivity. 


a  system  detection  efficiency  ?7spd  ~  70%  at  a  wave¬ 
length  of  1550  nm  with  dark  count  rates  of  roughly  4  Hz. 
A  noise-equivalent  intracavity  phonon  occupancy  ("Unep) 
can  be  defined  by  dividing  the  noise  count  rates  by  the 
per-phonon  sideband  photon  count  rate,  which  yields 


^nep  = 
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where  g  is  the  total  efficiency  of  the  setup,  including  the 
system  efficiency  of  the  SPDs  as  well  as  optical  insertion 
loss  along  the  path  from  cavity  to  detectors,  and  A  is  the 
transmission  of  the  filters  at  the  pump  frequency  relative 
to  the  peak  transmission.  The  above  equation  makes 
clear  the  benefits  of  large  cavity-enhanced  optomechani¬ 
cal  coupling  go ,  both  in  terms  of  the  low  power  sensitivity 
limited  by  detector  dark  counts  and  the  high  power  sen¬ 
sitivity  limited  by  pump  bleed-through.  Further  details 
about  the  detection  setup,  optical  and  mechanical  spec¬ 
troscopy,  and  derivation  of  tt-nep  can  be  found  in  the 
App.  [A]|F]  below. 

In  Fig.  it  we  show  the  measured  sensitivity  of  the 
phonon  counting  setup  for  A  =  cjm  (solid  grey  circles),  as 
well  as  the  expected  theoretical  curve  given  by  Eqn.[lJ  To 
measure  the  total  noise  count  rate,  the  pump  beam  is  de¬ 
tuned  far  off-resonance  from  the  optical  cavity,  A  u;m, 


FIG.  3.  a,  Phonon  count  rate  (blue)  and  inferred  phonon 
occupancy  (n)  (red)  as  a  function  of  intracavity  photon  num¬ 
ber  for  A  =  —  ujm-  Dashed  lines  indicate  points  below  (blue), 
at  (green),  and  above  (magenta)  threshold,  b,  Noise  power 
spectral  densities  (NPSD)  corresponding  to  the  dashed  lines 
in  a.  c,  Phase  plots  of  the  in-phase  (/)  and  in-quadrature  (Q) 
amplitudes  of  the  detected  optical  noise  corresponding  to  the 
dashed  lines  in  a 


such  that  negligible  motional  sideband  photons  are  gen¬ 
erated.  The  total  71nep  can  then  be  computed  by  nor¬ 
malizing  the  measured  count  rate  at  A  >  cjm  by  the 
per-phonon  sideband  count  rate  at  A  =  uom.  Here  the 
per-phonon  sideband  count  rate  is  extrapolated  from  the 
measured  sideband  count  rate  at  low  nc,  where  back- 
action  is  negligible  and  (n)  is  equal  to  the  room  temper¬ 
ature  thermal  mode  occupancy  of  no  ~  1100.  The  mea¬ 
sured  sensitivity  follows  the  expected  curve  at  low  power 
due  to  detector  dark  counts  (solid  red  curve),  but  at  high 
nc  saturates  to  a  value  several  times  larger  than  expected 
for  the  filter  suppression  of  the  pump  (solid  blue  curve). 
In  order  to  better  understand  this  excess  noise,  Fig.  [2p 
shows  measurements  of  the  tt-nep  as  a  function  of  filter- 
pump  detuning,  Amter?  a^  a  high  power  where  the  pump 
transmission  dominates  the  noise  (nc  ~  65).  A  strong  de¬ 
pendence  on  A fqter  is  observed,  with  a  peak  in  the  noise 
at  5  GHz  and  a  secondary  peak  at  6.1  GHz,  consistent 
with  laser  phase- noise  of  our  pump  laser  [26].  With  the 
addition  of  a  C-band  bandpass  filter  prior  to  the  SPD  to 
remove  broadband  spontaneous  emission  from  the  pump 
laser,  and  at  frequencies  far  from  the  laser  phase- noise 
peaks,  the  measured  unep  agrees  well  with  the  theo¬ 
retical  predictions  based  on  the  filter  pump  suppression 
(horizontal  dashed  curve).  At  the  relevant  detuning  of 
Amter  =  (vertical  dashed  curve),  we  measure  a  limit¬ 
ing  sensitivity  of  unep  =  0.89  db  0.05. 

In  what  follows,  we  focus  on  measurements  made  with 
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a  blue-detuned  pump  (A  =  —  u;m),  in  which  the  op¬ 
tomechanical  back-action  results  in  instability  and  self¬ 
oscillation  of  the  acoustic  resonator  0  [26] .  The  Stokes 
sideband  count  rate  detected  on  a  single  SPD,  shown 
versus  nc  in  Fig.  UK  displays  a  pronounced  threshold, 
with  an  exponential  increase  in  output  power  beginning 
at  nc  «  1200,  where  C  =  |7oM|/7i  ~  0.8,  in  agree¬ 
ment  with  the  expected  onset  of  instability  around  C  m  1 
(7  =  0).  This  sharp  oscillation  threshold  can  also  be  ob¬ 
served  from  the  measured  NPSD  (Fig.  [3)3) ,  in  which  the 
amplitude  of  the  mechanical  spectrum  is  seen  to  rapidly 
increase  with  a  simultaneous  reduction  in  linewidth,  and 
in  plots  of  the  in-phase  and  in-quadrature  components  of 
the  photocurrent  fluctuations,  which  show  a  transition 
from  thermal  noise  to  a  large  amplitude  sinusoidal  oscil¬ 
lation.  Also  shown  in  Fig.  S*  is  the  inferred  phonon  oc¬ 
cupancy  (n),  determined  by  calculating  the  per-phonon 
sideband  scattering  rate  at  low  nc,  where  (n)  is  roughly 
equal  to  the  thermal  occupation,  and  extrapolating  to 
higher  powers  assuming  linearity  in  the  optomechani¬ 
cal  interaction.  As  detailed  in  App.  0  a  self-consistent 
determination  of  the  oscillation  amplitude,  based  upon 
comparison  of  the  sideband  count  rates  above  and  below 
threshold,  indicates  that  even  at  our  highest  pump  power 
the  mechanical  amplitude  remains  small  enough  that  the 
linear  approximation  remains  valid. 

The  statistical  properties  of  the  resonator  near  the  self¬ 
oscillation  threshold  can  also  be  characterized  by  measur¬ 
ing  photon  correlations  using  a  HBT  setup  as  shown  in 
Fig-13  As  the  oscillation  threshold  is  crossed,  the  state  of 
the  acoustic  resonator  will  transition  from  a  thermal  state 
into  a  displaced  thermal  state  (DTS),  and  the  normalized 
phonon  intensity  correlation  function  near  r  =  0  should 
show  a  transition  from  bunching  (g^2\0)  >  1)  to  Pois- 
sonian  statistics  (g^2\r)  =  1  for  all  r).  Plots  of  g^2\r) 
below,  at,  and  above  threshold  are  shown  in  Fig.ffli.  Be¬ 
low  threshold  bunching  is  clearly  visible,  with  =  2 

as  expected  for  a  purely  thermal  state.  In  Fig.  Up  g(2)(0) 
is  plotted  versus  nc  for  both  blue-  and  red-detuned  pump 
light.  For  a  blue-detuned  pump  a  smooth  decrease  from 
s(2)( 0)  =  2  to  5(2)(0)  =  1  is  observed  in  the  threshold 
region,  while  for  a  red-detuned  pump  the  oscillator  is  ob¬ 
served  to  remain  in  a  thermal  state  through  threshold 
and  beyond.  The  decay  rate  of  the  acoustic  resonator, 
measured  from  both  the  linewidth  of  the  NPSD  and  from 
an  exponential  fit  to  g^2\r)  below  threshold,  is  plotted 
in  Fig.  |§-  The  decay  rate  as  measured  from  the  NPSD, 
which  includes  both  phase  and  amplitude  fluctuations, 
is  seen  to  increase  around  threshold  before  continuing  to 
decrease.  This  behavior  is  commonly  observed  in  semi¬ 
conductor  lasers  where  a  coupling  exists  between  the  gain 
and  the  cavity  refractive  index,  and  a  similar  effect  arises 
in  optomechanical  oscillators  due  to  the  optical  spring  ef¬ 
fect  m-  The  decay  rate  measured  from  g^2\r):  on  the 
other  hand,  which  measures  intensity  fluctuations,  be¬ 
gins  to  deviate  from  the  measured  linewidth  in  the  vicin¬ 
ity  of  threshold.  Thermal  phonon  emission  dictates  a 
strict  correspondence  between  the  second-order  and  first- 
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FIG.  4.  a,  Normalized  second-order  intensity  correlation  func¬ 
tion  g^2\r)  for  A  =  — cjm,  shown  below,  at  and  above  thresh¬ 
old  (left  to  right,  respectively).  Green  lines  show  a  simple  ex¬ 
ponential  fit,  while  black  lines  indicate  the  expected  theoret¬ 
ical  curve  using  decay  rates  measured  from  fitting  the  NPSD 
linewidth.  b,  Phonon  correlation  at  zero  time  delay  versus  nc 
for  A  —  —  (jjm  (blue)  and  A  —  ujm  (red).  The  top  and  bottom 
dashed  lines  indicate  the  expected  values  for  purely  thermal 
or  displaced  thermal  states  (DTS),  respectively,  c,  Mechani¬ 
cal  decay  rate  versus  nc  for  A  =  — o;m,  determined  from  the 
measured  linewidth  of  the  NPSD  (circles)  and  from  the  expo¬ 
nential  fit  to  g^2\r)  (diamonds),  d,  Fano  factor  versus  nc. 


order  coherence  functions  [4],  however,  above  threshold 
where  the  phonon  statistics  are  no  longer  purely  ther¬ 
mal,  such  a  deviation  is  possible,  and  in  fact  predicted 
for  self-sustaining  oscillators  m-  The  Fano  factor,  de¬ 
fined  as  F  =  (A n)2/(n)  =  1  +  (n)(g(2)(0)  —  1),  provides 
additional  statistical  information  about  the  fluctuations 
of  the  oscillator,  and  is  useful  for  defining  a  precise  os¬ 
cillator  threshold  [32]  as  well  as  distinguishing  between 
states  that  may  have  similar  or  identical  values  of  g ^  (0) 
(e.g.  a  coherent  state  versus  a  DTS)  [30].  The  Fano 
factor  of  our  mechanical  oscillator,  computed  from  the 
measured  g(2)(0)  and  the  inferred  values  of  (n),  is  dis¬ 
played  in  Fig.  [4]T  and  shows  the  expected  increase  and 
peak  in  fluctuations  at  threshold.  Above  threshold,  the 
Fano  factor  drops  again  due  to  saturation  in  the  optome- 
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chanical  gain,  approaching  a  measured  value  consistent 
with  that  expected  for  a  DTS  (F  ~  2nb  +  1). 

Although  we  have  emphasized  the  analogy  between 
the  optomechanical  oscillator  studied  here  and  a  laser, 
there  are  unique  differences  which  arise  due  to  the  in¬ 
trinsically  nonlinear  nature  of  the  radiation  pressure  in¬ 
teraction  in  an  optomechanical  cavity.  Recent  theoretical 
studies  Eol  1301133]  indicate  that  a  laser-driven  optome¬ 
chanical  oscillator  will  enter  a  nonclassical  mechanical 
state  with  anti-bunched  phonon  statistics  (F  <  1),  and 
under  slightly  more  restrictive  conditions,  strongly  neg¬ 
ative  Wigner  density.  Surprisingly,  this  is  predicted  to 
be  observable  even  for  classical  parameters,  i.e.,  outside 
the  single-photon  strong-coupling  regime  (go/ft  <  1),  and 
in  the  presence  of  thermal  noise.  Beyond  phonon  cor¬ 
relation  spectroscopy  of  optomechanical  oscillators,  it  is 
envisioned  that  sensitive  photon  counting  of  the  filtered 
motional  sidebands  may  be  utilized  in  the  preparation 
and  heralding  of  non- Gaussian  quantum  states  of  a  me¬ 
chanical  resonator  [6] .  For  the  OMC  cavities  of  this  work, 
with  their  large  optomechanical  coupling  rate  and  near 
millisecond-long  thermal  decoherence  time  at  sub-Kelvin 
temperatures  [28],  the  phonon  addition  and  subtraction 
processes  of  Ref.  [6j  should  be  realizable  with  high  fidelity 
and  at  rates  approaching  a  megahertz.  Whether  for  stud¬ 


ies  of  the  quantum  behavior  of  mesoscopic  mechanical 
objects  or  in  the  context  of  proposed  quantum  informa¬ 
tion  processing  architectures  utilizing  phonons  and  pho¬ 
tons  01,  such  photon  counting  methods  are  an  attrac¬ 
tive  way  of  introducing  a  quantum  nonlinearity  into  the 
cavity-optomechanical  system. 
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FIG.  5.  FEM  simulations,  a,  Electric  field  Ey  component 
of  the  optical  mode  at  frequency  cjc/27r  =  194  THz  (polariza¬ 
tion  in  the  plane  of  the  page  and  transverse  to  the  long  axis 
of  the  nanobeam),  b  Displacement  field  of  the  mechanical 
breathing  mode  at  frequency  cjm/27r  =  5.1  GHz. 
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Appendix  A:  Optical  and  mechanical  design 

The  optomechanical  crystal  (OMC)  studied  in  this 
work  is  numerically  optimized  for  both  optomechanical 
coupling  and  optical/acoustic  quality  factor,  via  finite- 
element  method  (FEM)  simulation  in  COMSOL  Mul¬ 
tiphysics  [34],  according  to  the  procedure  outlined  in 
Ref.  [25].  The  holes  on  the  ends  of  the  beam  support  si¬ 
multaneous  bandgaps  for  optical  light  with  wavelengths 
near  1550  nm  and  acoustic  waves  with  frequencies  from 
3  —  5.5  GHz,  while  the  variation  of  holes  towards  the 
center  of  the  beam  perturb  the  bandgaps  so  as  to  create 
co- localized  optical  and  acoustic  midgap  resonances.  The 
fundamental  optical  mode  (Fig.  HD  has  a  nominal  wave¬ 
length  of  1545  nm  and  the  fundamental  acoustic  mode 
(Fig.  [5Jd)  has  a  nominal  resonance  frequency  of  5.1  GHz. 
The  nominal  optomechanical  vacuum  coupling  rate,  due 
predominantly  to  photoelastic  effects,  is  predicted  to  be 
go/2ir  =  860  kHz.  Physically,  this  coupling  rate  is  the 
optical  resonance  frequency  shift  due  to  the  zero-point 
fluctuations  of  the  acoustic  resonator. 


Appendix  B:  Fabrication 

The  devices  are  fabricated  from  a  silicon-on-insulator 
(SOI)  wafer  (SOITEC,  220  nm  device  layer,  3/im  buried 
oxide)  using  electron  beam  lithography  followed  by  re¬ 
active  ion  etching  (RIE/ICP).  The  Si  device  layer  is 
then  masked  using  a  ProTEK  PSB  photoresist  to  de¬ 
fine  a  mesa  region  of  the  chip  to  which  a  tapered  lensed 
fiber  can  access.  Outside  of  the  protected  mesa  region, 
the  buried  oxide  is  removed  with  a  plasma  etch  and  a 
trench  is  formed  in  the  underlying  silicon  substrate  us¬ 
ing  tetramethylammonium  hydroxide  (TMAH).  The  de¬ 
vices  are  then  released  in  hydrofluoric  acid  (49%  aqueous 
HF  solution)  and  cleaned  in  a  piranha  solution  (3-to-l 
H2S04:H202)  before  a  final  hydrogen  termination  in  di¬ 
luted  HF.  In  fabrication,  arrays  of  the  nominal  design 
shown  in  Fig.  [5]  are  scaled  by  ±2%  to  account  for  fre¬ 
quency  shifts  due  to  fabrication  imperfections  and  disor¬ 
der. 


Appendix  C:  Experimental  setup 

The  full  experimental  setup  for  phonon  counting  and 
intensity  interferometry  is  shown  in  Fig.  [6]  A  fiber- 
coupled,  wavelength  tunable  external  cavity  diode  laser 
is  used  as  the  light  source,  with  a  small  portion  of  the 
laser  output  sent  to  a  wavemeter  (A-meter)  for  frequency 
stabilization.  The  remaining  laser  power  is  sent  through 
an  electro-optic  phase  modulator  (0-m),  used  to  gener¬ 
ate  optical  sidebands  for  locking  the  filter  cavities,  and 
a  variable  optical  attenuator  (VO A)  to  allow  control  of 
the  input  power  to  the  cavity.  The  signal  is  then  sent 
into  an  optical  circulator  which  sends  the  optical  probe 
to  the  a  lensed  fiber  tip  for  end-fire  coupling  to  the  de¬ 
vice.  The  cavity  reflection  can  then  be  switched  to  one 
of  two  detection  setups.  The  first  allows  the  signal  to 
be  switched  to  a  power  meter  (PM)  for  measuring  the 
reflected  signal  power  or  to  an  erbium-doped  fiber  am¬ 
plifier  (EDFA)  followed  by  a  high-speed  photodetector 
(PD).  The  resulting  photocurrent  can  be  sent  to  a  real¬ 
time  spectrum  analyzer  (RSA)  in  order  to  measure  the 
noise  power  spectral  density  (NPSD)  of  the  optical  signal 
or  to  a  vector  network  analyzer  (VNA)  which  can  be  used 
to  measure  the  full  complex  response  of  the  optical  cav¬ 
ity.  The  second  detection  setup  sends  the  cavity  reflec¬ 
tion  through  a  series  of  narrowband  tunable  Fabry-Perot 
filters  (~  50  MHz  bandwidth,  ~  20  GHz  free-spectral 
range)  in  order  to  reject  the  pump  frequency.  The  sig¬ 
nal  then  travels  through  a  variable  coupler  (VC)  and  is 
sent  to  the  dilution  refrigerator  where  is  is  detected  by 
two  superconducting  single  photon  detectors  (SPD).  The 
output  of  these  detectors  is  sent  to  a  time-correlated  sin¬ 
gle  photon  counting  (TCSPC)  module  for  calculation  of 
the  detection  correlation  function. 

Since  the  pump  laser  is  tuned  to  a  motional  side¬ 
band  during  the  phonon  counting  measurement,  the  two 
Fabry-Perot  filters  used  in  this  work  must  be  tuned  to 
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FIG.  6.  Experimental  setup.  Phonon  counting  setup.  A- 
meter:  wavemeter,  FPC:  fiber  polarization  controller,  0-m: 
electro-optic  phase  modulator,  VOA:  variable  optical  atten¬ 
uator,  SW:  optical  switch,  PM:  optical  power  meter,  EDFA: 
erbium-doped  fiber  amplifier,  PD:  fast  photodiode,  RSA:  real¬ 
time  spectrum  analyzer,  VNA:  vector  network  analyzer,  VC: 
variable  coupler,  SPD:  superconducting  single  photon  detec¬ 
tor,  TCSPC:  time-correlated  single  photon  counting. 


the  optical  cavity  resonance  via  an  initial  lock  and  sta¬ 
bilization  procedure.  This  procedure  is  as  follows.  First, 
the  pump  wavelength  is  tuned  to  the  blue  or  red  OMC 
sideband  by  optimizing  the  mechanical  transduction  sig¬ 
nal  on  a  spectrum  analyzer.  Since  the  power  of  the  ra¬ 
diated  Stokes-  or  anti- Stokes-scattered  light  is  too  low 
to  provide  a  feedback  signal  for  filter  stabilization,  we 
then  bypass  the  OMC  and  phase-modulate  the  pump  at 
frequency  cjm.  The  filters  are  then  locked  to  maximize 
transmission  of  the  sideband  which  is  resonant  with  the 
cavity.  After  a  stabilization  period  of  a  few  seconds,  the 
filter  positions  are  held  without  further  feedback  while 
the  pump  modulation  is  turned  off,  pump  power  is  ad¬ 
justed  for  the  desired  nc,  and  the  OMC  is  switched  back 
into  the  optical  path.  Once  locked,  the  transmission  of 
the  filters  is  observed  to  be  stable  to  within  5  —  10%  for 
several  minutes  in  the  absence  of  active  feedback  locking. 

In  order  to  avoid  pile-up  artifacts  [35]  in  the  acquired 
g (r)  histograms  presented  in  the  main  text,  the  photon 
count  rate  incident  upon  the  SPDs  is  kept  at  or  below 
30  kHz.  This  is  accomplished  with  a  variable  optical 
attenuator  on  the  output  of  the  filters,  and  is  sufficient  to 
maintain  a  flat  histogram  background  over  a  5  /is  window. 
The  absolute  count  rate  reported  in  Fig.  takes  the 
variable  attenuation  into  account. 


Appendix  D:  Device  characterization 

Full  characterization  of  the  optical  resonance  involves 
determination  of  the  single  pass  fiber- to- waveguide  cou¬ 
pling  efficiency  rjcp\,  the  total  energy  decay  rate  ft,  and  the 
waveguide- cavity  coupling  efficiency  r\K  =  fte/ft  (fte  is  the 
decay  rate  into  the  detection  channel).  The  fiber  collec- 


FIG.  7.  Optical  cavity  response.  a,  Optical  reflec¬ 
tion  spectrum  of  the  cavity  resonance  versus  cavity  detun¬ 
ing  A  =  ujc  —  u)\  (blue)  with  a  Lorentzian  fit  (red)  yielding 
Q c  =  236,  000.  b,  Phase  response  of  the  optical  resonance, 
yielding  fte/ft  =  0.52. 


tion  efficiency  is  determined  by  observing  the  calibrated 
reflection  level  far-off  resonance  with  the  cavity  and  is 
found  to  be  77cpi  =  0.63.  The  total  cavity  decay  rate  is 
determined  by  fitting  the  optical  reflection  spectrum  of 
the  cavity  (Fig.  [T^i)  and  yields  k,/2tt  =  818  MHz  (opti¬ 
cal  quality  factor  Qc  =  236,000).  The  reflection  level  on 
resonance,  when  normalized  to  the  off-resonance  reflec¬ 
tion  level,  is  related  to  the  cavity- waveguide  coupling  effi¬ 
ciency  by  R0  =  (1  —  2tjk)2.  However,  for  single-sided  cou¬ 
pling  this  is  not  a  single- valued  function  of  the  coupling 
efficiency.  Consequently,  the  complex  response  of  the 
cavity  is  measured  by  locking  the  laser  off-resonance  from 
the  cavity  and  using  a  vector  network  analyzer  (VNA)  to 
drive  an  electro-optic  modulator  (EOM)  and  sweep  an 
optical  sideband  across  the  cavity.  By  detecting  the  re¬ 
flected  power  on  a  high-speed  photodetector  connected 
to  the  VNA  input,  the  phase  response  of  the  cavity  can 
be  measured  (Fig.  Fh)  Fitting  this  with  prior  knowledge 
of  the  cavity  resonance  frequency  and  decay  rate  yields 
rjK  sat  0.52. 

To  characterize  the  acoustic  resonance,  the  cavity  re¬ 
flection  is  sent  through  an  erbium-doped  fiber  ampli¬ 
fier  (EDFA)  and  detected  on  a  high-speed  photodetec¬ 
tor.  The  EDFA  is  used  to  amplify  the  signal  so  that 
the  optical  noise  floor  overcomes  the  detector’s  electronic 
noise,  and  the  noise  power  spectral  density  (NPSD)  of 
the  optical  reflection  is  measured  on  a  real-time  spec¬ 
trum  analyzer  (RSA),  where  a  Lorentzian  response  due 
to  transduction  of  the  acoustic  thermal  Brownian  mo¬ 
tion  can  be  observed  at  the  acoustic  resonant  frequency 
cjm/27r  =  5.6  GHz.  For  a  pump  laser  locked  onto  the  red 
or  blue  mechanical  sideband  of  the  cavity  (A  =  uoc  —  uj\  = 
±cjm)  the  linewidth  of  the  transduced  signal  is  given 
by  7  =  7;  ±  7qm,  where  7om  =  ±kglnc/n  (rac  is  the 
steady  state  intracavity  photon  number).  The  depen¬ 
dence  of  linewidth  for  both  detunings  versus  nc  is  shown 
in  Fig.  [8j  By  averaging  the  two  sets  of  data  we  can  ex¬ 
tract  the  intrinsic  acoustic  damping  rate  7i/27r  =  3  MHz 
(Qm  =  1850),  which  is  seen  to  remain  constant  as  a 
function  of  nc.  By  fitting  the  excess  optomechanically 
induced  damping  70M  as  a  function  of  nc  we  extract  a 
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FIG.  8.  Calibration  of  go .  Mechanical  linewidth  7  versus  in¬ 
tracavity  photon  number  nc  for  A  —  ujm  (red)  and  A  —  —  ujm 
(blue).  The  intrinsic  linewidth  of  the  acoustic  resonator  71 
(black)  is  determined  by  averaging  the  blue  detuned  data  and 
yields  Qm  =  1850.  The  inset  shows  the  optomechanically 
induced  damping  70M,  obtained  by  subtracting  71  from  7, 
versus  nc.  The  linear  fit  shown  in  red  yields  a  vacuum  op¬ 
tomechanical  coupling  rate  of  go  =  645  kHz. 


coupling  rate  of  go  =  645  kHz. 


Appendix  E:  Single  photon  detectors 

The  detectors  used  in  this  work  are  amorphous  WSi- 
based  superconducting  nanowire  single-photon  detectors 
(SNSPDs,  or  SPDs  hereafter)  developed  in  collaboration 
between  the  Jet  Propulsion  Laboratory  and  NIST.  The 
SPDs  are  designed  for  high-efficiency  detection  of  individ¬ 
ual  photons  in  the  wavelength  range  A  =  1520  —  1610  nm 
with  maximum  count  rates  of  about  25  x  106  counts  per 
second  (c.p.s).  (reset  time  tR  =  40  ns)  [29]  and  very  low 
dark  count  rates  (DCRs).  With  the  SPD  in  the  supercon¬ 
ducting  state  (below  its  critical  temperature  Tc  =  3.7  K), 
a  DC  bias  current  I b  of  a  few  microamps  is  maintained 
through  the  nanowire  by  an  external  current  source.  The 
operating  range  of  the  detector  lies  between  a  lower  cutoff 
current  and  an  upper  switching  current  Jsw,  above  which 
the  nanowire  switches  to  a  non-superconducting  state. 
The  choice  of  quiescent  operating  current  if*  for  each  de¬ 
tector  was  made  to  roughly  maximize  the  ratio  of  the 
quantum  efficiency  ?7spd  to  DCR  while  operating  within 
the  ”  plateau”  region  where  both  77 spd  and  the  DCR  are 
nearly  constant  as  Jb  is  adjusted  slightly. 

The  SPDs  are  mounted  on  the  still  stage  of  a  3He/4He 
dilution  refrigerator  at  700  mK.  Single-mode  optical 
fibers  (Corning  SMF-28)  are  passed  into  the  refrigerator 
through  vacuum  feed-throughs  and  coupled  to  the  SPDs 
via  a  fiber  sleeve  attached  to  each  SPD  mount.  Proper 


alignment  of  the  incoming  fiber  with  the  15  /im  x  15  /am 
square  area  of  the  SPD  nanowire  is  ensured  by  a  self- 
aligned  mounting  system  incorporated  into  the  design  of 
the  SPD  [29  .  The  radio- frequency  output  of  each  SPD  is 
amplified  by  a  cold-amplifier  mounted  on  the  50  K  stage 
of  the  refrigerator  as  well  as  a  room-temperature  ampli¬ 
fier,  then  read  out  by  a  triggered  PicoQuant  PicoHarp 
300  time-correlated  single  photon  counting  module.  The 
counting  module  is  triggered  by  input  pulses  reaching  a 
voltage  above  a  fixed  discriminator  value  VJ-  Amplified 
photon-detection  pulse  heights  of  150  —  250  mV  are  typ¬ 
ical,  and  corresponding  discriminator  values  in  the  range 
110  —  150  mV  were  chosen  for  each  SPD  by  measuring 
nominal  count  rates  as  a  function  of  VJ  and  choosing  an 
operating  value  of  Vd  to  be  near  the  center  of  the  plateau 
region  in  which  the  observed  count  rates  are  independent 
of  small  changes  in  the  discriminator  setting. 

Initial  characterization  of  the  SPDs  was  centered  on 
measuring  dark  count  rates  and  the  quantum  efficiencies 
of  the  detectors.  The  measured  DCRs  are  sensitive  to 
various  channels  by  which  stray  light  may  couple  into 
the  fiber-detector  system,  including  ambient  laboratory 
lighting  and  thermal  radiation  both  inside  and  outside 
the  refrigerator.  By  tightly  spooling  (~  1.5  inch  diame¬ 
ter)  the  optical  fiber  within  the  fridge  to  filter  out  long 
wavelength  blackbody  radiation  and  systematically  iso¬ 
lating  the  optical  fiber  from  environmental  light  sources 
we  have  achieved  DCRs  of  2  —  4  c.p.s. 

Quantum  efficiency  measurements  of  the  SPDs  were 
made  using  laser  light  of  A  =  1554  nm  attenuated  to 
an  input  power  of  1.53  fW  at  the  input  to  the  fridge, 
corresponding  to  N  «  12,  000  incoming  photons  per  sec¬ 
ond.  We  calculate  ?7spd  by  referring  the  detected  photon 
count  rate  (less  the  corresponding  known  DCR)  to  the 
nominal  input  flux  N.  This  efficiency  incorporates  the 
intrinsic  detection  efficiency  of  the  SPDs  as  well  as  any 
losses  in  the  fiber  run  within  the  fridge  and  in  the  cou¬ 
pling  between  the  fiber  and  the  SPD  itself.  At  just  below 
the  respective  switching  currents  of  the  detectors,  we  find 
?7spd  =  70%,  with  this  result  depending  on  photon  po¬ 
larization  (<  20%  variability). 


Appendix  F:  Phonon  counting  sensitivity 

For  sufficiently  weak  optomechanical  coupling  (gQ 
k)  and  small  mechanical  amplitude,  the  equations  of 
motion  for  the  optomechanical  system  can  be  linearized 
about  a  large  steady  state  optical  field  amplitude.  For  a 
sideband  resolved  system  (k/2  cjm)  and  a  red-detuned 
pump  (A  «  cjm)  the  output  optical  field  may  then  be 
written  in  the  Fourier  domain  (in  a  frame  rotating  at  the 
pump  frequency)  as  [26] 
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where  ain(a;)  =  aS(cj)  -\-aVSLC(cj)  (a  is  the  steady-state  op¬ 
tical  field  at  the  pump  frequency,  avac(cj)  is  the  vacuum 
noise  of  the  pump) ,  /q  =  ft  —  fte  is  the  intrinsic  loss  rate 
of  the  optical  cavity,  d[(u:)  is  additional  vacuum  noise  ad¬ 
mitted  via  the  intrinsic  loss  channels,  and  b(u)  is  the  an¬ 
nihilation  operator  for  the  acoustic  resonator.  Note  that 
P(cj)  takes  the  place  of  b(cj)  for  a  blue-detuned  pump 
(A  «  -Um). 

As  6(uj)  is  sharply  peaked  around  c o  =  cjm,  we  can  spec¬ 
trally  filter  out  the  strong  optical  pump  at  uj  =  0.  The 
additional  optical  noise,  assumed  to  be  white  Gaussian 
noise,  cannot  be  filtered  out  in  this  way.  However,  in  the 
case  that  the  optical  noise  is  pure  vacuum  noise  it  will 
not  contribute  to  any  photon  counting  events.  Thus,  for 
the  purposes  of  photon  counting  the  output  optical  field 
can  be  written  post-filtering  as 


flout (t)  ~  2^Kencg°  b(t)  =  A  Av7om (F2) 
ft  V  ft 

which  shows  explicitly  that  in  this  linearized  regime 
photon  counting  is  equivalent  to  phonon  counting 
((aLt«out)  oc  (tfb)). 

As  can  be  seen  in  the  above  equation,  the  optically  in¬ 
duced  acoustic  damping  rate  7om  physically  represents 
the  per-phonon  rate  of  sideband  photon  emission,  corre¬ 
sponding  to  phonon  absorption  (emission)  for  A  =  co>m 
(A  =  —  o;m).  Of  the  sideband  photons  emitted  into  the 
optical  cavity,  a  fraction  fte/ft  are  subsequently  emitted 
into  the  detection  channel  and  detected  with  overall  sys¬ 
tem  efficiency  7,  including  both  the  system  efficiency  of 
the  SPDs  as  well  as  insertion  loss  along  the  path  from 
cavity  to  detector.  The  count  rate  per  phonon  is  thus 
given  by  Tsb,o  =  7(fte/ft)7oM5  and  the  total  count  rate 
is  given  by 


Ttot  =  rsB,o  (n)  +  Tpump  T  Tdark,  (F3) 


Tpump  .  Tdark 

^NEP  =  ^pump  +  ^dark  =  "W - b  - - •  (F4) 

J-  SB,0  1  SB,0 

The  dark  count  rate  Tdark  is  simply  a  measured  con¬ 
stant,  while  rpump  =  rjANpximp,  where  A  is  the  transmis¬ 
sion  of  the  filters  at  the  pump  frequency  relative  to  the 
peak  transmission,  and  7Vpump  is  the  input  photon  flux 
of  the  pump,  which  is  nearly  perfectly  reflected  from  the 
cavity  when  the  pump  is  far  off-resonance.  For  a  pump 
detuning  from  the  cavity  of  A  =  cdm,  the  input  photon 
flux  can  be  related  to  the  intracavity  photon  number  nc 
by  Apump  ~  ^mnc/^e-  Thus,  we  can  write  the  total 
noise-equivalent  phonon  number  as 
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Appendix  G:  Oscillation  amplitude  above  threshold 

The  classical  equation  of  motion  for  the  optical  cav¬ 
ity  field  a  is  given  in  the  frame  rotating  at  the  pump 
frequency  by 


b  —  —  (A  +  gox )  +  —  ^  a  +  H,  (Gl) 

where  x  is  the  position  of  the  acoustic  resonator,  H  = 
yf KeP[n/hujc,  and  Pin  is  the  optical  input  power  at  the 
device.  In  the  regime  of  parametrically  driven  self¬ 
oscillation,  the  amplitude  of  the  acoustic  oscillator  can 
become  large  enough  that  the  usual  linearization  ap¬ 
proximation  becomes  invalid.  However,  using  the  ansatz 
that  the  mechanical  oscillation  amplitude  is  given  by 
x(t)  =  /isin(c<;m£)  (note  that  (3  is  given  in  units  of  the 
zero-point  amplitude,  so  that  /32  =  4(n)  +  2),  an  exact 
solution  for  the  optical  cavity  field  can  be  written  as  a 
sum  of  sidebands. 

In  particular,  the  equation  of  motion  can  be  formally 
integrated  to  yield 


(G2) 


where  z  =  goP/u^.  This  can  be  solved  exactly  by  using 
the  Jacobi-Anger  expansion  etzc°8(d)  =  ^2ninJn(z)eine , 
where  Jn  is  a  Bessel  function  of  the  first  kind,  to  yield 


where  (n)  is  the  average  phonon  occupancy  of  the  acous¬ 
tic  resonator,  rpump  is  the  count  rate  due  to  residual 
pump  transmission  through  the  filters,  and  Tdark  is  the 
intrinsic  dark  count  rate  of  the  SPD. 

To  assess  the  sensitivity  of  the  phonon  counting  mea¬ 
surement,  the  noise  count  rate  can  be  divided  by  the  per- 
phonon  sideband  count  rate  to  obtain  a  noise-equivalent 
phonon  number 
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which  can  be  reindexed  to  have  the  form  a(t)  = 
aneinuJmt,  with 
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an  =  inn  Y,  Jm(*)  Jm-n(^,  (G4) 

m  m 

where  hm  =  k/2  +  z(A  +  rnu)m).  Note  that  if  we  are  only 
interested  in  the  total  energy  in  the  cavity,  |<a(£)|2,  the 
final  expansion  of  the  global  phase  factor  in  Eqn.  |G3|  is 
unnecessary  and  a  simpler  form  for  an  involving  no  sums 
may  be  used  [30,  36  .  However,  for  the  case  considered 
here  where  a  specific  frequency  component  of  the  cavity 
output  is  filtered  before  detection  it  is  necessary  to  use 
the  full  expression  given  in  Eqn.  |G4|  In  the  regime  of 
self-oscillation,  the  oscillation  amplitude  (3  is  determined 
by  balancing  the  optically  induced  mechanical  gain  with 
the  intrinsic  mechanical  loss 


7om 


Jn(z)Jn+l(z) 
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The  presence  of  the  filters  at  the  cavity  output  en¬ 
sures  that  the  SPD  count  rate  is  proportional  to  the 
number  of  intracavity  photons  in  the  first  Stokes  side¬ 
band  at  uj  =  u)\  —  o;m,  which  is  given  by  ni  =  |oq|2. 
For  the  lowest  input  power,  below  threshold,  n i  is  given 
by  the  simple  linear  approximation  rt\  =  [(4G2/ft)rib]/ft, 
where  rib  ~  1100  is  the  thermal  occupancy  of  the  me¬ 
chanical  resonator.  Near  and  above  threshold,  ni  can  be 
determined  from  the  ratio  of  the  above  and  below  thresh¬ 
old  SPD  count  rates  and  the  known  value  of  ni  below 


threshold.  For  a  given  input  power  with  n\  constrained, 
Eqs.  G4  and  G5 ,  along  with  the  condition  that  |7om  |  ~  7i 
above  threshold,  can  be  used  to  solve  for  the  detuning  A 
and  amplitude  (3  commensurate  with  the  observed  count 
rate. 


For  our  highest  input  power,  we  find  A  —  1.067cjm 
and  2  «  0.15.  This  amplitude  is  small  enough  that  the 
linear  approximation  (aq  oc  z)  is  still  valid.  In  particu¬ 
lar,  in  the  linear  regime  the  relation  Ji(z)  ~  zj 2  should 
hold.  For  the  largest  value  of  z  in  our  measurements,  we 
find  that  J\(z)  differs  from  z/2  by  only  about  0.3%.  The 
shift  in  detuning,  and  the  concomitant  reduction  in  oscil¬ 
lation  amplitude,  is  expected  due  to  the  thermo-optic  ef¬ 
fect  which  will  tend  to  shift  the  cavity  resonance  to  lower 
frequencies  as  the  total  intracavity  photon  population  is 
increased  by  the  amplified  Stokes  scattering  m- 

In  the  linearized  regime,  nc  (the  intra-cavity  pho¬ 
ton  number  of  the  0th  sideband  at  uq)  is  to  a  good 
approximation  equal  to  the  intra-cavity  photon  num¬ 
ber  in  the  absence  of  optomechanical  coupling,  nc  = 
(4ftePin/fi<^i)/(ft2  +  4A2).  The  total  optomechanical 
back-action  rate  (7om)  is  also  approximately  equal  to 
the  scattering  rate  from  the  0th  sideband  to  the  1st 
Stokes  sideband  (70,1)  in  the  linearized  regime,  which  for 
A  =  —  uom  blue  detuned  pumping  yields  |7om|  ~  70,1  ~ 
4 G2 /k.  As  this  is  the  case,  in  the  main  text  we  don’t  dif¬ 
ferentiate  between  back-action  damping  and  Stokes  scat¬ 
tering  rate,  and  nc  is  a  simple  placeholder  for  Pin  in  all 
of  our  plots. 


